The typical mass ratio and typical final spin in supermassive black hole mergers 
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We prove that merging supermassive black holes (SMBHs) typically have neither equal masses, 
nor is their mass ratio too extreme. The majority of such mergers fall into the mass ratio range 
of 1 : 30 to 1 : 3, implying a spin flip during the inspiral. We also present a simple expression for 
the final spin Xf of the emerging SMBH, as function of the mass ratio, initial spin magnitudes, and 
orientation of the spins with respect to the orbital plane and each other. This formula approximates 
well more cumbersome expressions obtained from the fit with numerical simulations. By integrating 
over all equally likely orientations for precessing mergers we determine a lower approximant to the 
final spin distribution as function of the mass ratio alone. By folding this with the derived mass 
ratio dependent merger rate we derive a lower bound to the typical final spin value after mergers. 
We repeat the procedure deriving an upper bound for the typical spin in the case when the spins are 
aligned to the orbital angular momentum, such that there is no precession in the system. Both slopes 
of Xf as function of the initial spins being smaller than one lead to two attractors at ^ rec = 0.2 and 
^aiign _ q 4^ respectively. Real mergers, biased toward partial alignment by interactions with the 
environment (accretion, host galaxy, etc.) would generate a typical final spin lying between these 
two limiting values. These are the typical values of the spin after the merger, starting from which 
the spin can built up by further gaseous accretion. 
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I. INTRODUCTION 

Einstein's general theory of relativity predicts that the 
coalescence of two compact objects (neutron stars or 
black holes) is accompanied (and driven by) intense grav- 
itational radiation. Stellar mass (a few to a few ten solar 
masses, Mq) black hole binaries emit gravitational waves 
with frequency falling into the best sensitivity range of 
LIGO [if, Virgo @ and GEO600 [3] Earth-based inter- 
ferometric gravitational wave detectors. Up to date there 
are very few observations Q indicating the existence of 
intermediate mass black holes. Binaries formed by such 
black holes would emit gravitational waves falling into the 
frequency range of third generation gravitational wave 
detectors, like the Einstein Telescope (4J. 

Supermassive black holes (SMBHs) with masses of 
10 6 -7- 3 x 10 9 M Q (or perhaps even higher) on the other 
hand are quite frequent, residing in the centre of each 
sufficiently massive galaxy. Their growth occurs by ac- 
cretion phases and by mergers, the estimated contribu- 
tion of each of these processes to the growth of mass 
being model-dependent. Recent observations @ show 
that some galaxies never merge, and yet may have cen- 
tral back holes; representing either the pure SMBH birth 
population, or the birth population with some gaseous 
accretion. 

The accretion process has been modeled with the in- 
clusion of magnetic fields, electromagnetic radiation of 



the disk and energetic jets transporting angular momen- 
tum from the polar regions [7|-[9|. Beside increasing the 
mass, accretion will also spin up the black holes. The 
spin limit reached due to accretion by canonical black 
holes (the system of a black hole and electromagnetically 
radiating accretion disk), as expressed in terms of the 
dimensionless spin is \can = 0.998, close to the theo- 
retically allowed maximum for a black hole, the unity. 
Indeed, such a high spin powering the jets seems com- 
pulsory for understanding the low energy cutoff in the 
energetic electron spectra of jets in radio galaxies [Toj . 
Active Galactic Nuclei, in particular the closest, Cen A 
are the most likely sources of for the Ultra High Energy 
Cosmic Rays [TT| . 

When galaxies merge, eventually their central SMBHs 
will also do so. Dynamical friction transfers some of the 
orbital angular momentum of the binary black hole sys- 
tem to the stellar environment, being ejected at the poles, 
a process which drives the system through the last par- 
sec [HI- Other mechanisms to overcome the last parsec 
are relaxation processes due to cloud/star - star interac- 
tions, repopulating the stellar orbits in the center of the 
galaxy [l3j , binary orbital decay by three-body interac- 
tions in the gravitationally bound stellar cusps [HI, or 
the interplay of three accretion disks: one around each 
black hole and the third, circumbinary, removing orbital 
angular momentum from the binary [15fl . 

At about 0.005 parsecs gravitational radiation takes 
over dynamical friction as the leading dissipative effect 
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(l6| . For many SMBH binaries the gravitational waves 
emitted in the process of coalescence fall into the fre- 
quency range the long-delayed space mission LISA [l7j . 
Depending on how rich in gas the binary environment 
may be and whether there any circumbinary disk has 
been formed, certain alignment between the proper spins 
and the orbital angular momentum could occur due to 
the Bardeen-Petterson effect (based in turn on the Lense- 
Thirring precession) [l8j . The two situations which could 
occur are the mergers precessing under random angle 
(also known as dry mergers) and non-precessing mergers, 
implying complete alignment of the spins and orbital an- 
gular momentum (wet mergers). The randomness in the 
orientation of precessing mergers typically reduces the 
final spin [19]. For equal mass precessing mergers this 
varies from 0.69 for non-spinning black holes up to val- 
ues of 0.73 for maximally spinning black holes [20] - For a 
mass ratio 1 : 10 the range of final spins opens up to the 
interval between 0.2 and 0.83, respectively, as function 
of the initial spins. For non-precessing mergers all con- 
figurations would practically conserve a high initial spin 
during the inspiral. 

When the two black holes are at large distances, the 
orbital angular momentum L is always much larger then 
the individual spins Si- However at the characteris- 
tic radius r* w 0.005 parsec 1 (and the corresponding 
post-Newtonian (PN) parameter e* = Gm/c 2 r* rj 10~ 3 , 
defining the beginning of the inspiral, see (HI), where 
gravitational radiation takes over dynamical friction as 
the leading order dissipative effect, the ratios Si/L ss 

(e*) 1 ^ 2 q 3 ~ 2l Xi depend on the mass ratio q > 1. For a 
maximally spinning larger black hole and separation r* 
the ratio is one at about q 30. For mass ratios larger 
than 30 therefore the spin dominates over the orbital an- 
gular momentum during the whole inspiral. 

During the inspiral gravitational radiation further re- 
duces the orbital angular momentum, but not the spin 
magnitudes. The spins will only precess driven by the 
leading order spin-orbit coupling and corrections due to 
spin- spin and mass quadrupole - mass monopole coupling 
{32J ■ In the process the direction of the total angular mo- 
mentum remains unchanged, in an averaged sense over 
one radial orbit [33| . 

Gravitational radiation does not modify this conclu- 
sion on short time-scales. Radiative evolutions with spin- 
orbit [3~3|, spin- spin [35j and mass quadrupole - mass 
monopole [36j couplings have been investigated, and their 
analysis in Ref. [37J lead to the important result that the 
instantaneous radiative changes of the spins average out 
during a radial period. On this timescale therefore there 



1 The distance r* depends weakly, as m 5 / 11 on the total mass and 
negligibly, as rj 2 / 11 on the symmetric mass ratio rj = [i/m = 
(q 1 / 2 + g-i/2) - 2 ! where 

fi = mim2/m is the reduced mass. 



is no secular radiative change of the spin vectors at all: 



This result confirms that the spin dynamics can be re- 
garded as a pure precession, up to high PN orders in- 
cluding radiation reaction, on the timescales comparable 
with a radial orbit. On much larger timescales however 
the spin vector will undergo a reorientation (spin- flip) , as 
explained in detail in Ref. [l6j |. 

There are several scenarios possible, according to the 
actual mass ratio:: 

a) The masses are comparable mi w m±. In this case 
at the end of the inspiral (when the PN approximation 
breaks down) the orbital angular momentum still domi- 
nates over the spins (l6| . [2l| . Radiating away this rem- 
nant orbital angular momentum during the merger phase, 
while extrapolating the conservation of the direction of 
the total angular momentum and of the individual spin 
magnitudes to this phase [HI , (3(| could significantly re- 
duce the final spin in all cases when the individual spins 
were severely misaligned with each other and with the 
orbital angular momentum. Such a misalignment would 
be typical in the case of precessing mergers. Therefore 
for equal masses in a precessing merger a not too high 
final spin can be considered typical. 

b) The mass ratio is in the range 1 : 30 to 1 : 3. In this 
case the orbital angular momentum dominates over the 
spin only at the beginning of the inspiral, and as such is 
roughly aligned with the total angular momentum. At 
the end of the inspiral however the orbital angular mo- 
mentum becomes smaller than the dominant spin, which 
has therefore to be reoriented towards the invariant to- 
tal angular momentum direction. For precessing mergers 
this process causes a spin-flip during the inspiral, but 
does not reduce significantly the magnitude of the domi- 
nant spin [l6| , [2ij ! Non-precessing mergers on the other 
hand already imply an alignment of the spins and orbital 
angular momentum, therefore neither the spin magni- 
tude, nor its direction will be changed by this process. 
In both cases, whatever happens to the orbital angular 
momentum during the plunge, its small value (compared 
to the dominant spin) at the end of the inspiral will ob- 
struct any serious further change in the final spin. 

In this mass ratio range therefore the magnitude of the 
dominant spin will not be much reduced by the merger, 
nevertheless a significant reorientation of its direction 
during the inspiral will typically occur for precessing 
mergers, which could be followed only by a minor fur- 
ther spin-flip during the plunge. 

c) The mass ratio is less than 1 : 30. Then the or- 
bital angular momentum is too small from the beginning 
of the inspiral to modify the dominant spin. Neither its 
magnitude, nor its direction are affected and we practi- 
cally face the inspiral of a test mass into the much larger 
black hole. 

In this paper we revisit the merger process, based on 
the recent data of Ref. [22|. We first derive the SMBH 
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mass distribution. In Section |TT] we fit a broken power 
law for the differential mass function, then, based on this 
fit and a number of simple and reasonable assumptions 
we derive the mass ratio distribution. We note that the 
results of Ref. [22| are fully consistent with earlier results 
based on much smaller statistics [icj . 

Next we derive in Section IIIII a simple approximant 
for the final spin of the emerging SMBH, as function of 
the mass ratio, initial spin magnitudes, and orientation 
of the spins with respect to the orbital plane. In the 
Appendix we compare the approximant with the more 
cumbersome expressions existing in the literature, which 
were obtained by fit to numerical simulations. 

In Subsection IIV Al we adopt the configuration of pro- 
cessing mergers, which allow for all relative spin and 
orbital angular momentum orientations on equal foot- 
ing, lowering the chances for a large final spin after the 
merger. By integrating over all orientations in the pro- 
cessing merger limit (without allowing any preference for 
alignment), for any initial spin set we determine a lower 
approximant to the final spin distribution as function of 
the mass ratio alone. By folding this with the previously 
derived mass ratio dependent merger rate, we obtain a 
lower bound to the typical final spin after SMBH mergers. 

By contrast, in the non-precessing merger limit there 
is a perfect alignment of the spins with the orbital angu- 
lar momentum, hence the integration should be carried 
on for this configuration alone, and only over the mass 
ratios, folded with the mass ratio distribution. By this 
method, in Subsection IIV Bl we get an upper bound for 
the typical final spin. 

We discuss the implications of our results and present 
the concluding remarks in Section IVl 



II. MASS RATIOS IN SMBH MERGERS 

In Ref. [lj] we gave a simple preliminary estimate of 
the typical mass ratio of merging SMBHs. We revisit the 
problem more rigorously here, both from a mathemat- 
ical point of view and by employing new, more precise 
data on SMBH masses, presented in Ref. (22[. We do 
note that selection effects strongly influence some statis- 
tical arguments, in the case, that selection is based on 
detectable activity at the center of a galaxy for instance, 
on a far-infrared or ultra-violet excess; in the first case 
this could be due to selecting for central emission lines, in 
the second due to a central star-burst, and in the third to 
a visible central accretion disk. Our approach, taken in 
this paper, does suffer from the selection effect, that the 
work done by [22| used the colors of an old stellar popu- 
lation as the starting point, and then cut the sample to 
include only early Hubble type galaxies. However, allow- 
ing for a sample of late Hubble type galaxies would not 
increase the merger rate very much, since such galaxies 
usually suffer few if any mergers [6[ . 



A. The differential mass function 
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FIG. 1: (Color online) The integral mass function represented 
on logarithmic scale shows a remarkably good fit to the data 
from Ref. The differential mass function is taken as a 

broken power law with powers —I and —3, the breakpoint 
being at 8.9 x IO 7 M . 



The SMBH distribution <S>bh{M B h) can be inter- 
preted as a power law with an exponential cutoff [22| . 
This can be well approximated by a broken power law 
[H-[23, also confirmed by the survey [H. The SMBH 
integral mass function data represented on Fig 5. of 
Ref. [22j . after omitting the first two data points which 
do not refer to black holes, but rather to nuclear star 
clusters [27|, also suggest the differential mass func- 
tion &bh{Mbh) oc Mgff, with a — 1, starting from 
the lower mass limit of m a w 10 6 M Q to the break- 
point, which is approximately at m* w 10 8 Mq; then 

^bh(^bh) oc Mg^, with f3 = 3, starting from m* to 



the upper mass limit, taken here as nib rj 3 x 10 M Q . 
We prove this statement in the remaining part of the 
subsection. 2 

The SMBH data is represented on Fig[TJ which shows 
the integral mass function 3> {Mbh) dMsH (in 

Mpc -3 ) as a function of the SMBH masses (in M©), rep- 
resented on log-log scale. Due to the breakpoint (and by 
normalizing all masses to m+), the integral mass function 



2 Note, that the lower mass data points were not important for the 
considerations in Ref. [2^| . concerned mainly with the highest 
energy cosmic rays, such that a different fit of acB = 2 was 
advanced there. Nevertheless the supermassive black holes with 
lower mass are important in the merger statistics, therefore in 
this paper we chose a = 1 due to the tendency of the first data 
points to be aligned horizontally (see FigfTJ. The limit m a is 
lowered here as compared to the choice of Ref. [HI such that the 
mass of the SMBH in the centre of our Galaxy is not the lower 
mass limit any more. 
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for any Mbh < m* is 
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while for any Mbh > m* is, respectively 
$ (M BH ) dM BH 
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Here A: is a dimensional normalization constant. Both 
expressions reduce to km*/ 2 at Mbh = m *- Comparing 
with the data at m* allows to fix log (km*/2) w —3, thus 
km*/ 2 = 1(T 3 . Therefore 



log 



<S>{M B H)dM B H 



Mb 



, 1 + 4.6 (log - x) . if M B h < m* ■ 

10g( lQ-2x m 2 ifM B H>m* ) ' (4) 



where x = log Mbh- Because Mbh is given is solar 
masses, so is m*. The broken power law with powers —1 
and —3 gives the best fit with the data by setting the 
breakpoint at m* — 10 7 - 95 M Q w 8.9 x 10 7 M Q , as seen on 
Fig [U (Note that the breakpoint turns out to be shifted 
as compared with the number given in Ref. [16].) The 
fit is remarkable, the sum of the squares of the deviances 
between the points and the function values, divided by 
the square of the error bars is only 0.22. 



B. SMBH mass ratio distribution 

Based on the new SMBH mass function derived in the 
previous subsection here we work out the estimates for 
the likelihood of the mass ratios, following the logic of 
Ref. [16j. However the changed mass values imply more 
cases to be included in the analysis. 

The number of encounters for a given mass ratio q = 
TO1/TO2 > 1, represented as dM / dq is proportional to the 
product of the distribution functions for both black holes, 
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FIG. 2: (Color online) The number of SMBH encounters with 
mass ratios q as function of log 2 q. 



folded with the merger probability F and integrated over 
the mass m2 of the smaller black hole:: 



dAfjg) 
dq 



oc 



m b /q 



$ , BH(rn 2 )§ BH(qm 2 )F(q, m 2 )dm 2 



(5) 

The merger probability in turn is proportional to the 
cross section (we neglect the weak dependence on the 
relative velocity of galaxies as these are not too high, the 
Universe being not old enough for mass segregation). 

In order to determine the cross section, we assume that 
each galaxy merger is followed by the merger of their 
central SMBHs. Therefore we basically evaluate the cross 
section of merging galaxies. Further, we note that the 
masses of the galaxies and their central SMBHs correlate 
due to 

• the correlation of the mass of the central SMBH 
with the mass of the host galactic bulge [281, 

• the proportionality of the mass of the central 
SMBH with both the spheroidal galaxy mass com- 
ponent and the total mass (including dark matter) 
of the galaxy [29j. 

It is likely that the more massive SMBH, thus the 
most massive galaxy dominates the cross section, thus 
the merger rate F. As the cross section is a function of 
the galaxy mass (thus SMBH mass), we take F ~ (qm 2 ) ■ 
We chose £ = 1/2 based on the following observation: 

• the comparison of our galaxy with dwarf 
spheroidals shows that an increase by a factor of 
10 in radius (thus 10 2 in cross section) is accom- 
panied by an increase by a factor of 10 4 in mass 



31]. 



The break point m* splits the SMBH range into two 
intervals, encompassing a mass range of about a factor 
of qi — 89 and q 2 = 36. Thus (by normalizing all masses 
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to m*) we estimate for any q £ [1,36] the number of 
encounters as 



dM(q) 
dq 



gg[l,36] 
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TO* 
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m 2 q 



-a / \ f 
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,*/g V m */ \ m * / V m * / 

!>/<? / tt;, 2 \ ^ / TO2<?\ ^ / ?Tl2(?\^ 



* / I m * J \ m * J 



dm 2 ■ (6) 



The first, second and third lines of the right hand side of 
Eq. © contain, respectively, the mergers of: two SMBHs 
from the lower mass interval; one SMBH from the lower, 
the other from the higher mass interval; and both SMBHs 
from the upper mass interval. The condition q < q 2 
assures that the upper limit of the integrals is larger than 
the lower limit. 

For q £ [36, 89] there are two contributions, arising 
from the combination of either two light or a light and a 
heavy SMBHs: 

dJV(q) 

—j — <x 

a 1 q£ [36,89] 



:n, \ ( TO2Q\ / m 2 q \ 

1* J \ m± ) 2 



m* 



IH ' m 2 



k / q \TO* J \ TO* J \ TO* J 

The condition q < qi again assures that the upper limit 
of the integrals is larger than the lower limit. 

Finally for q £ [89, 3000] there is one single contribu- 
tion 

dJV(q) 



d<l ,je[89,3000] 



m b /q 



m 2 
m* 



/ m 2 q \ I m 2 q \ 



to* 



V m* 



dm 2 . (8) 



Eq. ([7]) expresses the encounters of a light SMBH from 
the lower interval with a heavy SMBH from the upper 
interval. 

Integration over m 2 gives 



<W(g) 
dq 



<3G[1,36] 



1 + C - 2<5 



1 + £. - 2/3 



(9) 



and 



dJV(q) 



dq 



q£ [89,3000] 



1+Z-a-P 



(11) 

(We have employed to^/to* = q 2 and to*/to q = qi.) 
With the preferred parameter values 5 = 1, @ = 3 and 
£ = 1/2, qi = 89, q 2 = 36 the number of encounters for 
the three ranges simplifies to 
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Here we have normalized such that 

36 AKfl ^ ,89 dM{g) 



(12) 



dq 



dN{q) 
dq 9 e[i,36 

3000 dN{q) 
>'■> dq ge [89,3000] 



dq 



dq g g [36,89] 

dq = 1 (13) 



holds. 

Defining the number of encounters in a mass ratio in- 
terval [q\,q 2 ] as N qi ^ q2 = ^-jq^-dq we obtain the 
percentages of the mergers with the mass ratio ranges 
falling between [1,3], [3,30], [30,100] and [100,3000], re- 
spectively as 



M+3 = 12-1 % 
Mo-ioo = 29.2%, 



A/3-30 = 48.9 % , 
A/ioo +3000 — J-o % . (14) 



The distribution of the mass ratios is shown in more de- 
tail on the histogram of Fig [2] 

The most likely mass ratio range, occurring in approx- 
imately half of the mergers, turns out to be q £_(3 , 30), 
in agreement with the rough estimate of Ref . [16| . This 
is the mass ratio, where a spin flip occurs during the in- 
spiral [l6j . The second most numerous mass ratio range, 
approximately in 30% of the cases is for q £ (30, 100). 
Both the comparable mass case with q £ (1,3) and the 
extremal mass ratio case, defined here as q £ (100, 3000) 
represent just about 10% each of the SMBH mergers. 
This important result makes compulsory to model SMBH 
mergers for non-equal masses. 



dM(q) q- 1+& - q^- i+2a q^- a 



dq 
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(10) 



III. AN APPROXIMATE FINAL SPIN 
FORMULA IN SMBH MERGERS 

In this section we propose a formula for the final spin, 
which on the one hand approximates reasonably well 
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more cumbersome expressions derived from fits with nu- 
merical runs, on the other hand is simple enough to fa- 
cilitate the numerical integrations we will carry on in 
the remaining part of the paper. In this section we use 
v = q^ 1 . 

In the system with the Newtonian orbital angular mo- 
mentum on the z-axis and the periastron on the x-axis 
the spins are 

G 

Si = — m/^^ 2l_3 Xi (sin K.j cos Ci, sin/tj sin^j, cos Kj) . 
c 

The magnitude of the total spin S = Si + S2 is found 
from 



[Si + 2S X ■ S 2 + S|] 

G 

— mfj, 

c i=l,2 



21 V2 



2i-3 



Xi 



2XlX2 cos 7 



1/2 



with 



cos 7 = cos Ki cos K2 + sin k\ sin/i2 cos ((2 — Ci) (15) 

while the orbital angular momentum (assuming circular 
orbits, thus Ln — [irv), to leading order can be written 
as 

G 



-mps 



-1/2 



(0,0,1) 



where e = Gm/c 2 r = v 2 /c 2 is the post-Newtonian pa- 
rameter. (This parameter increases as the black holes 
approach each other.) 

The dimensionless version Z = cJ/ Gm/i of the magni- 
tude of the total angular momentum reads 

3 = 7^— [(L N + S).(L N + S)] 1/2 
= [L 2 V +2L N -(S 1 +S 2 ) + S 2 ] 1/2 



[e- 1 + 2e- 1 / 2 ]T i/^xicos* 

i=l,2 

+ E (^ 2 ^ 3 X J ) 2 + 2xiX 2 cos 7 ] 1/2 



The final spin magnitude is denoted 



Sf = 



G 



l fXf 



We identify as an upper limit for the final spin the 
magnitude of the total angular momentum at the end 
of the inspiral, obtaining \f — V (m/rrif) 2 J /. Here 
3/ = 3 (e = £/) and 77 = fi/m = v (1 + v) 2 . By in- 
troducing the efficiency of mass conversion into gravita- 
tional radiation as eaw = 1 — mf/m, we can express 
m/rrif = (1 — eaw) ■ Hence 
V_ 

Xf , 2 

(1 - ecw) i=1 2 



V 21 3 Xi COS Ki 



i=l,2 



Xi 



■ 2xi X2 cos 7] 



1/2 



(16) 



We are interested in establishing a lower boundary for 
the value of the final spin, thus we set the efficiency to 
zero. The maximal value of the bracket for any given \i 
arises when the spins are aligned with the orbital angular 
momentum: 

x max = ^ = v ^-1/2 + „-l Xi + ^ 



With e 



-1/2 



z.f — 2 (at two Schwarzschild radii, this is 
the radius of the innermost bound circular orbit in the 
Schwarzschild geometry) and for maximal spins this gives 
X™ = 1, irrespective of the actual value of v. Therefore 

we normalize Xf by setting e ^ =2 in Eq. ([T6|) . and 



obtain a very simple expression for the final spin: 
(1 + v? L 



Xf 



+ E(- 2 



V 2% 3 Xi COS Ki 



■ 2xiX2 cos 7] 



1/2 



(18) 



When the mass ratio is extreme (y — > 0), Eq. p^|) cor- 
rectly reproduces Xf = Xi> a result to be expected from 
the test particle limit. This final spin function qualita- 
tively reproduces well the more cumbersome final spin 
expressions found in the literature from fits with numer- 
ical runs. In the Appendix we compare in detail the ex- 
pression (fT5|) with the one presented in Ref . [39[ , finding 
that for the largest part of the parameter space it slightly 
underestimates the final spin. 



IV. THE TYPICAL FINAL SPIN 
A. Precessing (randomly oriented) mergers 

In this subsection we discuss the typical spin in the 
merger of two black holes by assuming generic precess- 
ing mergers. As this implies complete randomness in 
the relative angular momenta orientations, we integrate 
the final spin formula (|18[) over all possible orientations. 
Then we weight this orientation independent, but still 
mass ratio dependent final spin with the probabilities for 
a given mass ratio (|12p derived earlier in Section [H] and 
integrate over the mass ratios, obtaining a typical final 
spin as function of initial spin magnitudes only. As the 
integration over the mass ratio implies to integrate over q 
(according to the method of evaluating the merger rate), 
we will rewrite v = q^ 1 in all expressions. 



1. Mass ratio dependent typical final spin 

We first integrate the expression of the final spin (|Tg|) 
over all spin directions. By adopting the precessing 
merger model, we allow for random spin orientations. 
The assumption of randomness sets cos Ki and Q , the co- 
sine of the spin polar angles and the spin azimuthal angles 
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as evenly distributed random variables. Instead of the in- 
dividual azimuthal angles, the combination 7 = £ 2 — £1 
appearing in Eq. (|18j) and representing the relative spin 
azimuthal angle will be randomized. 

Integrating over all orientations (and properly normal- 
izing by 8n) we find therefore a lower bound for the mass 
ratio dependent final spin as: 

X P / eC (Xi,q) = -^J J J Xf (Xi,Q, K i,l)dcosK 1 dcosK 2 d'f 

(19) 

The final spin ^ rec for precessing mergers (random 
configurations) as function of xi — Xi an d log q is rep- 



resented in the left panel of Fig. |3l as arising from a 
numerical integration (with approximation by the mid- 
point method). For equal masses the final spin ranges 
from 0.5 (for nonspinning black holes) to 0.6 (for maxi- 
mally spinning black holes). This is consistent with the 
corresponding results of Ref. [lo| . For q > 100 we have 
the test particle limit: the final spin is accurately approx- 
imated by xi ( m other words the orbital angular momen- 
tum does not modify the spin of the larger SMBH). In 
between there is the mass range with the most frequent 
encounters. For q 10 for example the lower bound for 
the final spin ranges from 0.18 (nonspinning mergers) to 
0.85 (maximally spinning mergers). 





FIG. 3: (Color online) The typical final spin in supermassive black hole mergers as function of xi = X2 and log q, represented 
for precessing mergers (averaged over random configurations) - left panel; and mergers with the spins and orbital angular 
momentum fully aligned - right panel. 



2. Typical final spin in precessing mergers 

We establish an overall typical final spin for precessing 
mergers by integrating Eq. (|19[) over all possible mass ra- 
tios, weighted with the mass ratio dependent probability 
of encounters given in Eqs. (|12[) : 

prec { \ f prec / \ dj\f {^q) 

Xf {Xi) = / Xf KXhl) — dq 

J Ji 1 dq 3£[li36] 

+ / Xf fe>?) — ; dq 

J 36 a 1 q£ [36,89] 

+ xT c ^,q)^- m 

J 89 a Q q<E [89,3000] 

The result of the numerical integration can be seen as 
the lower curve in Fig 01 

We note that for merging SMBHs in fast rotation xi = 
X2 ~ 0.998 (the canonical spin limit, which occurs when 



both the accretion and the radiation of the disk are taken 
into account @) the final spin becomes Xf ~ 0.75. 



B. Non-precessing (aligned) mergers 

There is no precession in the perfectly aligned con- 
figurations, when the two spins and the orbital angular 
momentum are parallel. Such configurations could arise 
due to accretion or by other mechanisms. We do not 
model such mechanisms here, just assume the alignment 
of the spins and orbital angular momenta of the two-body 
system. 



1 . The mass ratio dependent final spin 

The spins being aligned to the orbital angular momen- 
tum implies k% = = 7. Inserting these values in Eq. 
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(I18p . the final spin for non-precessing mergers as a func- 
tion of the initial spin magnitudes and mass ratio takes 
a remarkably simple form: 



mass ratios 1 : 3 to 1 : 30, the interesting mass ratio 
range where a spin-flip would occur during the inspiral 
161. 



Xf (Xi, Q) = - 772 ( 2 + IXi + q- X X2) ■ (21) 
(l + q- 1 ) 

The final spin Xf %9n f° r non-precessing mergers is repre- 
sented in the right panel of Fig [3] as function of \i = X2 
and log q. 



2. Typical final spin in non-precessing mergers 

Next, we again integrate over the mass ratios, by prop- 
erly weighting with the mass ratio dependent probability 
of encounters, as given in Eqs. (1T21) : 



dq 



Xf (Xi) = / Xf (Xi,q) —, — 
1 Jl 1 dq qe[136] 



dM(q) 



J 35 a Q qe [36,89] 

^3000 



dq 



dAf{q) 



>m "9 q£ [89,3000] 

The result of the numerical integration can be seen as 
the upper curve in Fig 2] 

For merging SMBHs in fast rotation xi = Xi ~ 0.998 
(the canonical spin limit) the final spin is Xf ~ 0.86, 
much higher than for precessing mergers, however still 
reduced essentially as compared to the initial spin values. 



V. CONCLUDING REMARKS 

In this paper we have studied the typical mass ratio 
and the typical final spin in a two-body system com- 
posed of supermassive black holes (SMBH), thus we did 
not consider the perturbations induced by either of the 
accretion disks, nearby stellar population, magnetic fields 
or jets. SMBHs reside in the center of each galaxy and 
following the frequent galaxy mergers they also merge. 
Various dissipative processes, like dynamical friction, ac- 
cretion and emitted gravitational radiation are respon- 
sible to their gradual approach and finally gravitational 
radiation is which drives them to coalescence through a 
sequence of inspiral, merger and ringdown. 

By starting from precise and new data on the SMBH 
mass distribution we derived both a differential and an 
integral mass function, shown on Fig. [TJ The differential 
mass function is a broken power law, with coefficients —1 
and —3, with the breakpoint at 8.9 x 10 7 M Q . 

Then, exploiting a number of simple and reasonable 
assumptions we derived the mass ratio dependent prob- 
ability of encounters of two such SMBHs, represented on 
Fig. [U This confirms our expectation that the most 
frequent (approximately half of the) encounters are for 



aligned mergers 

precessing mergers 




00 0.1 02 03 0.4 05 0.6 0.7 08 09 10 



Ai — X2 



FIG. 4: (Color online) The typical final spin \t as func- 
tion of xi = X2 only, in the randomly precessing and the 
non-precessing merger limits (lower and upper curves, re- 
spectively). The curves are obtained by integration over all 
mass ratios of the expressions (|19p ^P rec (%i = X2> <?) and (|21[) 
xj'" s " (xi = X2,q), respectively, weighted with the mass ratio 



dependent probabilities of encounter (|12[1 . The line of equal 
initial and final spins is also indicated. Where this line crosses 
the final spin curves, there are two attractors (denoted by 
large dots), to where the final spin would converge after a 
sequence of mergers in the two scenarios. 

Next, based on certain well-founded assumptions we 
derived a simple analytical expression for the final spin 
of such a merger, depending on the mass ratio, initial spin 
magnitudes, and orientation of the spins with respect to 
the orbital plane and each other. This formula approxi- 
mates well more cumbersome expressions obtained from 
the fit with numerical simulations, and it is much simpler, 
thus advantageous in order to carry on the cumbersome 
numerical integrations which followed. 

We proceeded to find the typical final spin in two lim- 
iting and highly idealized scenarios. First we allowed 
for perfectly random orientations (precessing case), over 
which we have integrated, obtaining a final spin still de- 
pending on the initial spin magnitudes and mass ratio. 
Then we folded with the derived mass ratio dependent 
merger rate, we integrated over the mass ratio, deriving 
a lower bound to the typical final spin value after merg- 
ers. 

Secondly we considered the non-precessing configura- 
tion, with all spins and the orbital angular momentum 
perfectly aligned. Folding the final spin for this partic- 
ular configuration again with the derived mass ratio de- 
pendent merger rate and integrating over the mass ratio 
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we obtained an upper bound for the typical spin. These 
are represented as function of the initial spin magnitudes 
(chosen to be equal 3 ) on Fig. |U A third curve, the line 
of equal initial and final spins is also indicated on the 
figure. The fact that both slopes of Xf as function of 
the initial spins are smaller than one, leads to important 
consequences. 

If we imagine a sequence of idealized (either randomly 
precessing or non-precessing) mergers, what happens is 
that low spins tend to increase by mergers while high 
spins decrease. There are in fact two attractors at 
^rec _ q 2 anc j ^A h a n _ 0.45, respectively, where the 
spins converge after a reasonable number of the two types 
of mergers. 

Real mergers, biased toward partial alignment by in- 
teractions with the environment (accretion, host galaxy, 
etc.) would generate a typical final spin lying between 
these two limiting values. Indeed, for example the galaxy 
group distribution around NGC383 4 looks like a spin- 
dle, with the spin of the central black hole in the domi- 
nant galaxy, NGC383, aligned with the long axis of the 
spindle. This shows a correlation between the dominant 
galaxy black hole spin and the distribution of the other 
galaxies, all with central black holes as well. It is to be 
expected that the distribution of galaxies in the environ- 
ment of a dominant galaxy is not random, but correlated, 
such that in a merger of a galaxy with the dominant 
galaxy the final spin, depending on the nature of the 
correlation, could fall anywhere between the two curves 
shown of Fig. [U 

After the merger episode gaseous accretion can start 
to increase the spin again. If gaseous accretion were 
strong, then the spin could become quite large in rela- 
tively short time. We propose to work out quantitatively 
such a model in a forthcoming work. 
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Appendix A: Comparison of the final spin formula 
with related results 

Various papers have presented empirical formulae for 
the final spin, with the functional form partially moti- 
vated by PN expressions and coefficients fitted to the 



3 The second dimensionless spin parameter X2 will anyhow have 
but a small impact on the result, as the ratio of the total spins 
Si scales with the mass ratio squared. 

4 This is also known as the 3C31 radio source, cf. the NASA 
Extragalactic Database. 



result of numerical runs. In one of the latest such works 
Barausse and Rezzola [3!| suggested a set of criteria I-V) 
such a formula should obey. We compare this formula 
with our Eq. (fTS]). 

The condition I) of Ref. [39[ implies no rest mass loss 
by gravitational radiation; our lower spin limit estimate 
with ecw = does the same. Condition III) assumes 
that the radiation loss in the last stage of the merger is 
along the direction of the total angular momentum, this 
pro per ty continuing to hold similarly as during the inspi- 
ral [33j. By identifying \f with Z we also assume that. 
Similarly, this identification assures the validity of condi- 
tion IV) , which translates to no change in Ki and 7 during 
the plunge. However we note that in the strict sense the 
spin-spin and quadrupole-monopole couplings of the PN 
dynamics will obstruct this assumption; therefore this as- 
sumption cannot be considered valid for any distance, as 
assumed in Ref. [3{|. Nevertheless these angular evolu- 
tions are negligible on the orbital timescale, thus during 
the plunge (over which we assume its validity and which 
lasts only from a fraction of an orbit to a few orbits) the 
condition can be regarded as accurate. Condition V. of 
Ref. [39[ implies that the initial spins should drop out 
completely from the final spin formula in the particular 
case of equal masses and equal, but opposed spins. Our 
Eq. (fTB|) can be specified for this configuration by insert- 
ing v = 1, X2 = Xii cos 7 = — 1 and cos ft 2 = — coski 
and it gives ^-*P ec co ™/ l 9- = e~^ 2 /4, which is also inde- 
pendent of the initial spins. Therefore condition V) also 
holds. 

Finally, condition II) assumes that there are three vec- 
tors with conserved length. These are the two spins (that 
we also assume) , and the vector Ln — J + Sj (what we 
do not). The latter condition in our notation implies 

const. = (e- 1/2 L N - (3 - rj^Xf) j) 

= e- 1 + {2-r,- 1 Xf) 2 

-26- 1 / 2 (3 - rT 1 */) L N • J . (Al) 

We can set the constant to ej 1 by evaluating the for- 
mula at Sf, where Xf = V-Sf- Note that £ _1 a r thus it 
changes with f, at Keplerian order. For generic r we have 
X/=const. and 3 changing significantly only on the ra- 
diation timescale (due to gravitational radiation). Over 
the orbital timescale the change in 3 by gravitational ra- 
diation backreaction is at 2.5PN orders, while over the 
precessional timescale is of 1PN order. The evolution of 

a — cos -1 ^Ln • Jj generates a 1PN change over the or- 
bital timescale [4jJ, therefore the leading order changes 
over the orbital timescale of the second and third terms 
in the second line on the right hand side of Eq. (|A1|) are 
of order e 5 / 2 and e 1 / 2 . Thus we conclude that condition 
II) concerning the constancy of the length of the vector 
Ln — J + S/ cannot be extended to arbitrary r, as in fact 
changes with e 1 / 2 . 

In the most generic case discussed in Ref. (39j . their 
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Eqs. (6) and (8) reproduce our Eq. (IT6)) . provided we 
replace their |[|, given by their Eq. (10) and rewritten in 
our notations as 

2 

|[| = 2^3-3.5171 ^- + 2.5763 — — 



0.4537 7T f-^ - 0.8904 



+ - 



(i+f) 2 



l + v 2 



(Xl COS Ki + V 2 X2 COS K 2 ) 



0.1229 



(Xi + ^ 4 X2 + 2^ 2 XiX2COS 7 ) (A2) 



— 1/2 

with e j, ' . In what follows, we compare the two values 



161) for the final spin, once computed by replacing e 



-i/2 



with 



given by Eq. (jA2|) . then with e / 1/2 = 2. For 
example in the equal mass v = 1, equal spin %i = %2 case, 
when the spins are opposed to each other, thus 7 = tt 
and K2 = 7r — Ki, the ratio xf R /Xf is identically 1.37, 
regardless of the values of xi — Xi an d «i, therefore x/ 
underestimates X^ fl - Various other configurations, all 
for equal dimensionless spins %2 = Xii are represented 
on Figs. [5] and [HI 




FIG. 5: (Color online) The final spin estimates \f (green surfaces) and Xf (magenta) as function of \i = X2 and v for perfect 
alignment of the spins with the orbital angular momentum (non-precessing case - upper row); and anti-aligned spins in the 
plane of motion (severe precession - lower row). Except a narrow parameter range with high mass ratio and high spin values in 
the upper row configuration, the estimated Xf ls smaller than xf R - The agreement increases with decreasing v\ for the aligned 
configuration (upper row) is better in the high spin regime (visible on the right panel), then for low spin (left). 
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1.0 



FIG. 6: (Color online) The final spin estimates Xf (green surfaces) and xf R (magenta) as function of xi = X2 an d Ki for mass 
ratios v = 1 (upper row), v — 0.1 (middle row) and v = 0.01 (lower row). The represented configuration has the smaller spin 
confined to the plane of motion and the larger spin lying in the plane span by the smaller spin and orbital angular momentum. 
The agreement increases with decreasing v (with Xf < Xf aTi large v); is better in the high spin regime (visible on the right 
panel), then for low spin (left); is also better for configurations with «i 6 [0, 7r/2] than for the severely misaligned configurations. 
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